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7 It is shown that a three-dimensional (3D) modified-kinetic Alfven waves 

a (m-KAWs) can propagate in the form of Alfvenic tornadoes characterized 

9 by plasma density whirls or magnetic flux ropes carrying orbital angular mo- 

10 mentum (OAM). By using the two fluid model, together with Ampere's law, 

11 we derive the wave equation for a 3D m-KAWs in a magnetoplasma with m e /rrii <C 

12 f3 <C 1, where m e (m,) is the electron (ion) mass, (5 = A7rn kB(T e +Ti) / 'Bq, 

o no the unperturbed plasma number density, ks the Boltzmann constant, T e (T e ) 

a the electron (ion) temperature, and B the strength of the ambient magnetic 

is field. The 3D m-KAW equation admits solutions in the form of a Laguerre- 

16 Gauss (LG) Alfvenic vortex beam or Alfvenic tornadoes with plasma den- 

17 sity whirls that support the dynamics of Alfven magnetic flux ropes. 
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1. Introduction 

is The discovery of the Alfven wave Alfven [1942] had significantly changed the landscape 

lo of collective phenomena in magnetized plasmas that frequently occur in interstellar media 

20 Shukla and Dawson [1984] and in large laboratory plasma devices [Lennen et al. , 1999; 

21 Kletzing et al. , 2010; Gekelman et al. ,2011]. In the Alfven wave, the restoring force 

22 comes from the pressure of the magnetic fields, and the ion mass provides the inertia. 

23 The dispersion relation for the low-frequency (in comparison with the ion gyrofrequency) 

24 Alfven wave is deduced from the magnetohydrodynamic (MHD) equations [Alfven , 1942; 

2 5 Cramer , 2001], which do not account for the wave dispersion effects. The dispersion 

26 [Stefant , 1970; Brodin and Stenflo , 1990; Morales et al. , 1994; Morales and Maggs , 

27 1997; Shukla and Stenflo , 1999; Dastgeer and Shukla , 2009; Damino et al. , 2009] to 

28 the Alfven wave comes from the finite frequency (u/uj C i) and the magnetic field-aligned 

29 electron inertial force in cold magnetoplasmas, as well as from the finite ion gyroradius 
so effect and the gradient of the electron pressure perturbation in a warm magnetoplasma, 

31 Due to non-ideal effects Brodin and Stenflo [1990]; Morales et al. [1994]; Morales and 

32 Maggs [1997]; Shukla and Stenflo [1999]; Damino et al. [2009], Alfven waves thus 

33 couple to fast and slow magnetoacoustic waves, the inertial and kinetic Alfven waves 

34 Stefant [1970]; Morales et al. [1994]; Morales and Maggs [1997]; Dastgeer and Shukla 

35 [2009]; Damino et al. [2009] in a uniform magnetoplasma. The magnetic field-aligned 

36 phase speed (ou/k z ) of the inertial (kinetic) Alfven wave is much larger (smaller) than the 

37 electron thermal speed. Accordingly, the inertial (kinetic) Alfven wave arises in a plasma 

38 with j3 <C m e /mi (m e /m,i <C j3 <C 1). The dispersive shear Alfven wave plays a significant 
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39 role in the auroral Chaston et al. [2008a, b], magnetospheric Sundkvist et al. [2005], 

40 solar Tripathi and Gekeklman [2010]; Salem et al. [2012], astrophysical and laboratory 

41 [Gekelman et al. ,2011, 2012] magnetoplasmas with regard to acceleration and heating of 

42 the plasma particles Shukla et al. [1998], reconnection of magnetic field lines, wave-wave 

43 and wave-particle interactions [Hasegawa and Chen , 1976; Carter et al. , 2006; Shukla 

44 et al. , 2007], and the formation of nonlinear Alfvenic structures [Sundkvist et al. , 2005; 

45 Stasiewicz et al. , 2003; Shukla et al. , 2012]. Specifically, it is likely that the dispersive 

46 Alfven waves power the auroral activities Chaston et al. [2008a, b] and the solar coronal 

47 heating [Tomczyk et al. , 2007; Jess et al. , 2009; Mcintosh et al. , 2011], as well as may be 

48 responsible for the formation of Alfvenic tornadoes Widemeyer-Bohm et al. [2012] which 

49 can transport magneto-convective energy from one region to another in solar plasmas. 

so In this Brief Report, it is shown that a 3D m-KAW can propagate as a twisted vortex 

si beam in a plasma with m e /mj <C /3 <C 1. The present work thus complements a recent 

52 investigation Shukla [2012] that focused on examining a 3D twisted inertial Alfven wave. 

53 Twisted m-KAWs can be related with magnetic flux ropes or Alfvenic tornadoes of differ- 

54 ent scalesizes in observational data from laboratory experiments Gekelman et al. [2012], 

55 and also from the solar plasmas Widemeyer-Bohm et al. [2012]; Salem et al. [2012]. 

2. Mathematical description 

se We consider a magnetized electron-ion plasma in the presence of low-frequency (in com- 

57 parison with the electron gyrofrequency u ce = eB /m e c, where e is the magnitude of the 

58 electron charge, Bq the strength of the external magnetic field zBq, m e the electron mass, 

59 and c the speed of light in vacuum, z the unit vector along the z-axis in a Cartesian coor- 
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60 dinate system) m-KAWs with the electric and magnetic fields E = — V0 — (l/c)zdA z /dt 

61 and Bj_ = V A z x z, respectively, where and A z are the scalar and magnetic field- aligned 

62 vector potentials, respectively In the m-KAW fields, the electron density perturbation is 

63 obtained from 



dt Aire dz 
64 where we have used the electron fluid velocity 



dn el + ^dV^ = Q? 



u e ~^zxV^ - x Vn ei + z-^—VlA z , (2) 

±>o &o n o Alien® 

65 which uses Ampere's law that relates A z and the magnetic field-aligned electron fluid 

66 velocity u ez . Here n e i(^C n ) is a small electron density perturbation in the equilibrium 

67 density n . 

68 Since the parallel phase speed of the m-KAWs is much smaller than the electron thermal 

69 speed Vr e = (kBTe/rrie) 1 / 2 , one can neglect the electron inertia and obtain from the 

70 magnetic field aligned electron momentum equation 



dd> edA z kBT e dn e \ , . 

az c dt n oz 

?! which dictates that the parallel electric force —n§eE z and dpi/dz are in balance, where 

72 the magnetic field-aligned electric field E z = — <90 — c~ x dA z fdt and pi = k,BT e n e i is the 

73 electron pressure perturbation. 

7 4 We can now eliminate A z from Eqs. (1) by using (3), obtaining 
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n e i - 



4ne dz 2 



c 2 d 2 V\<P 



= 0, 



(4) 



where \dc = (ksTe/ Annoe 2 ) 1 / 2 is the electron Debye radius. 

The perpendicular (to z) component of the ion fluid velocity Ujj_ is determined from 



which is obtained by manipulating the perpendicular component of the ion momentum 
equation, where u C i = eB /miC is the ion gyrofrequency and the ion mass. We have 
assumed that the m-KAW phase speed is much larger than the ion thermal speed, and 
therefore neglected in Eq. (5) the contribution of the ion pressure gradient. The ions 
are assumed to be confined in a two-dimensional (x — y) plane perpendicular to z. The 
magnetic field-aligned ion dynamics is unimportant for the m-KAW s, since the phase 
speed of the latter is much larger than the ion-sound speed. 
From the linearized ion continuity equation and Eq. (5), we readily obtain 



where 7iji(<S n ) is a small ion number density perturbation. 

Invoking the quasi-neutrality condition n el = rin — ni, we eliminate rin from Eq. (6) by 
using Eq. (5), obtaining the wave equation for a 3D m-KAW after elementary calculation 






(6) 



d 2 



d 2 V 2 



dz 2 



C\ ( d 2 2 \ d 2 ] 



(7) 
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as where Ca = B / \fAnnomi is the Alfven speed, and p s = C s /u C i is the sound gyroradius, 

89 with C s = (k B T e /mi) l l 2 being the ion-sound speed. The effect of the ion temperature can 

90 be incorporated by re-defining p s = (C s /u C i)(l + 3Tj/T e ) 1//2 . 

91 Within the framework of a plane-wave approximation, assuming that n\ is proportional 

92 to exp(— iuut + ik ■ r), where u and k(= kj_ + zk z ) are the angular frequency and the wave 

93 vector, respectively, we Fourier analyze (8) to obtain the frequency spectra Shukla and 

94 Stenflo [1999] of the modified (by the u/u ci effect) KAWs 



2 k 2 z V 2 {\ + k\p 2 s] (R) 



95 where kj_ and k z are the components of k across and along z, and Aj = c/u pi the ion 

96 inertial scale length, with u pi = (47rn e 2 /mj) 1 / 2 being the ion plasma frequency. The 

97 /c 2 Aj-term in Eq. (8) comes from the u/u C i effect. 

98 Let us now study the property of a twisted m-KAW. We seek a solution of Eq. (7) in 

99 the form 



rii = N(r) exp(ik z z — lu^t), (9) 

100 where N{r) is a slowly varying function of z. Here r = (x 2 +y 2 ) l l 2 and k is the propagation 

101 wave number along the axial (the z— axis) direction. By using Eq. (10) we can write Eq. 

102 (7) in a paraxial approximation (viz. \d 2 N/dz 2 \ <^ k 2 N) as 



21 + v N = °' (10) 
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103 where Uk = k z VA/(l + k^X 2 ) 1 ^ 2 is the angular frequency of the magnetic field-aligned 

104 KAW. Here Z — (1 + k 2 \ 2 )k z z and is in unit of p s . Furthermore, we have denoted the 
us operator V^iV = (l/r)(d/dr)(rdN/dr) + (l/r 2 )d 2 N/d9 2 , and introduced the cylindrical 
wo coordinates with r = (r, 9, z). 

107 The solution of Eq. (10) can be written as a superposition of Laguerre-Gaussian (LG) 

108 modes [Allen et al. , 1992, 2003; Mendonca et al. , 2009], each of them representing a 

109 state of orbital angular momentum, characterized by the quantum number /, such that 



N = '£N p iF pl (r,Z)exp(iW), (11) 
P i 

no where the mode structure function is F p i(r, z) = H p iX^L^(X)exp(—X/2), with X = 
in R 2 /W 2 (Z), r = r/p s , and W = 2tt(1 + 47r 2 X 2 /L 2 z )z/L z , L z = 2-n/k z , and W(Z)/L Z 
n2 representing the normalized width of a twisted Alfvenic vortex beam. The normalization 
us factor H p i and the associated Laguerre polynomial L^(X) are, respectively, 



1/2 

(12) 
n4 and 



H p i = 



2V^ 



(1 + pV- 



L[\X) 



exp(X) d p 
Xp\ dXP 



X l+P exp(-X) 



(13) 



us where p and I are the radial and angular mode numbers of the DSAW orbital angular 
U6 momentum state. In a special case with / = and p = 0, we have a Gaussian beam. 
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U7 The LG solutions, given by Eq. (11), describe the structure of a twisted 3D Alfvenic 

us vortex beam carrying OAM. In a twisted m-KAW vortex beam, the wavefront rotates 

no around the beam's propagation direction (viz. the z-axis) in a spiral that looks like fusilli 

120 pasta (or a bit like a DNA double helix), creating a vortex and leading to the m-KAW 

m vortex (KAWV) beam with zero intensity at its center. A twisted m-KAWV beam (or 

122 an Alfven tornado) can be created with the help of two oppositely propagating 3D m- 

123 KAWs that are colliding in a magnetoplasma. Twisting of 3D m-KAWs occurs because 

124 different sections of the wavefront bounce off different steps, introducing a delay between 

125 the reflection of neighboring sections and, therefore, causing the wavefront to be twisted 

126 due to an entanglement of the wavefronts. Thus, due to angular symmetry, Noether 

127 theorem guards OAM conservation for a m-KAWV beam that is accompanies the parallel 

128 electric field, sheared magnetic field, and finite density perturbations. 

3. Summary and discussions 

129 In summary, we have shown that 3D m-KAWs in a uniform magnetoplasma can prop- 

130 agate in the form of Alfvenic tornadoes which have axial and radial extents of the order 

131 of the ion inertial length c/u pi and the ion sound radius p s , respectively. The Alfvenic 

132 swirls can be identified as observational signatures of rapidly rotating magnetic field ropes 

133 [Tripathi and Gekeklman , 2010; Gekelman et al. , 2012] (or magnetic field tornadoes), 

134 which can provide an alternative mechanism for channelling electromagnetic energy and 

135 heat fluxes from the surface of the Sun into the corona Widemeyer-Bohm et al. [2012] 
we through Alfvenic tornadoes/whirls. Furthermore, the present investigation of twisted m- 
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137 KAW vortex beams can also be useful for diagnostic purposes, since the m-KAW beam 

us frequency can be a fraction of the ion gyrofrequency. 
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